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§ : 1 Introduction 

O ■ The variance-covariance matrix of a multivariate normal distribution is usually estimated 
by the sample variance-covariance matrix, which is distributed as Wishart distribution. 
^ ! Let S be distributed according to Wishart distribution VFp(z/, S), where p (> 2) is the 
^ \ dimension, v (> p) is the degree of freedom, and S is the variance-covariance matrix of 
the original multivariate normal distribution. 

In many situations of multivariate analysis, such as principle component analysis, 
canonical correlation analysis, we need to estimate the eigenvalues of S rather than S 
itself. Also, many test statistics in multivariate analysis have distributions determined 
solely by the eigenvalues of E because of their invariance property under some natural 
transformations. 

For the estimation of the eigenvalues of S, the corresponding sample eigenvalues of 
S are usually used, but their distribution is quite complicated and makes it difficult 
to obtain mathematically clear results. Especially in a decision theoretic approach we 
encounter difficulty since we essentially need the calculation of the risk (the expectation 
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of a loss) with respect to the distribution of the eigenvalues for finite degrees of freedom 
V. Mainly because of this difficulty, there exist only a few literature which directly deal 
with the estimation of the eigenvahies from the standpoint of the decision theory. Dey 
(1988) and Jin (1993) derive estimators which dominate the traditional estimators under 
the (non-scale-invariant) quadratic loss function. In view of the decision theory, one of 
the important tasks is to derive an admissible estimator, but it has been an unsolved 
problem so far. The aim of this paper is the derivation of an admissible estimator. For 
the proof of admissibility, we adopted the method of Gosh and Singh (1968), in which 
they proved the admissibility of an estimator for the reciprocal of the scale parameter of 
Gamma distributions using "Karlin's method" (Karlin (1958)). 

Here we formally state the framework. Let Ai > . . . > Ap > denote the eigenvalues 
of S, while Zi > . . . > > are the eigenvalues of S. As is well known, the distribution 
of Z = (Zi, . . . , Zp) depends only on A = (Ai, . . . , Ap). For an estimator 

V>(0 = (^i(0,---,^pW), 

we measure the loss by the scale-invariant squared error loss function 

Y,m) - \f/>^i = Y.m)/>^^ - ^f- (i) 

i=l 1=1 



2 Main Result 

Before stating the main result as a theorem, we introduce some notation. For a vector 
X — (xi, . . . , Xp) and a set of powers a. — (cti, . . . , a^) the monomial denoted 
by x^. \l a = ai = ■ ■ ■ = Up is common, we denote the monomial by a;". Let H = 
(hij) denote a p-dimensional orthogonal matrix. The group of p dimensional orthogonal 
matrices is denoted by 0{p) and ii is the invariant probability measure on 0{p). Since 
we mainly work with the reciprocal of the population eigenvalue, tj = Xj^ {j = 1, . . . ,p), 
more often than A^ itself, we define the following notation for convenience. 

= {t = (ti,...,tp) I (0 <) a <ti < ■■■ <tp < 6 (< oo)}, 
G{l) = l(^-p-^y^l[{l,-l^), 

i<j 



F{t\l) = / exp ( -^EE^^^^-^i I M^)^ 



The density function f{l\t) oi I is given by 

mt) ^Kr/'G(l)F{t\l), (2) 
where X is a constant (not depending on I and t). Our main result is given as follows. 
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Theorem 1 For I < i < p, let 



The estimator = {ipl{l) , . . . , ip*{l)) is admissible with respect to the loss function 

Proof of this theorem is given in Section HI 
Notice that (1 < i < p) can be rewritten as 

^/^Ki) = j2^r,m, (3) 

where 

^ ^^''^ ^ U) {-\ ELi iskK,) tj dt,iH)dt ^ > 

It is easily seen that Tjj(Z) (1 <i,j < p) is scale-invariant, that is, Tij{cl) = Tij{l) for any 
positive constant c. Furthermore r^j's are nonnegative and 



Z/ + 2' 



since Ylj ^Ij ~ ^- This means that — (V^KO? • • • 5 "^^(0) is an estimator which 

shrinks + 2). 

i^iil) (1 < < p) has another useful expression; 



V 7=1 ^ 



(5) 



where 

fij{l) is again bounded and scale- invariant. (Lemma [H] in Section HI) 



3 Some simulation studies 

In this section, we report a small-scale simulation result which illustrates the behavior 
of the admissible estimator i/'*(Z) = (i'Kl), ■ ■ ■ ,ippil)) compared to the simple estimator 
cp*il) = l/{u + 2). 
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Table 1: p = 



2 





u = 


5 


u = 


20 


A 


risk of 


risk of (f* 


risk of i/)* 


risk of if* 


(1.0,1.0) 


0.623 


0.785 


0.186 


0.266 


(1.0,0.9) 


0.605 


0.756 


0.170 


0.230 


(1.0,0.8) 


0.557 


0.678 


0.150 


0.193 


(1.0,0.7) 


0.550 


0.636 


0.161 


0.188 


(1.0,0.6) 


0.565 


0.653 


0.169 


0.183 


(1.0,0.5) 


0.571 


0.634 


0.186 


0.182 


(1.0,0.4) 


0.559 


0.606 


0.189 


0.175 


(1.0,0.3) 


0.605 


0.643 


0.202 


0.183 


(1.0,0.2) 


0.591 


0.614 


0.189 


0.180 


(1.0,0.1) 


0.608 


0.627 


0.192 


0.189 



Using the variable transformation, 

ti = riu, t2 = r2U, . . . , tp_i = rp_iM, tp = u 
we can easily notice that (jl]) equals 

m ELi ELi rshhl,r^^^/^^^K^'-' r;^\f d^^iH)dr 

where r = (ri, . . . , rp_i) and 9^q = {r | < ri < ■ • ■ < rp_i < 1}. For given p, u, I, we 
calculated Tij{l) (1 < i,j < p) using 100 random points uniformly distributed respectively 
on 0{p) and 

We carried out a simulation for p = 2 and p = 3. In each case, z/ equals 5 or 20 and 
several patterns of population eigenvalues A are given. We used 1000 Wishart random 
matrices for the risk calculation for each p, i/, A. Table [T] and [2] show the simulation 
results. We notice that if the population eigenvalues are close to one another, then the 
estimator -0* performs better than c^*, while as population eigenvalues get dispersed, 
the risk of i/'* rapidly increases. This indicates the admissibility of i/j* is acquired by 
the good performance when population eigenvalues are equal at the expense of the poor 
performance when they are scattered. 

4 Proofs 

In this section we give a proof of Theorem [TJ Since the proof is long and complicated, 
first we give an outline of the proof for readability. Then we give a full proof in a series 
of lemmas. Long proofs of some lemmas are given in Appendix. 
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Table 2: p = 3 





u = 


5 


u = 


20 


A 


risk of "0* 


risk of ifi* 


risk of ■0* 


risk of (f* 


(1.0,1.0,1.0) 


0.860 


1.472 


0.254 


0.525 


(1,1/2,1/4) 


0.730 


1.071 


0.326 


0.286 


(1,1/3,1/9) 


0.859 


1.018 


0.726 


0.291 


(1.0,1.0,0.5) 


0.766 


1.286 


0.243 


0.363 


(1.0,0.5,0.5) 


0.761 


1.076 


0.263 


0.337 


(1.0,1.0,0.2) 


0.933 


1.244 


0.527 


0.377 


(1.0,0.2,0.2) 


0.979 


1.101 


0.348 


0.355 


(1.0,1.0,0.1) 


2.073 


1.213 


1.480 


0.362 


(1.0,0.1,0.1) 


1.838 


1.148 


0.756 


0.354 



An outline of the proof. Assume that some estimator '0(0 — (^1(^5 • • • ; V'p(O) dom- 
inates 0*(Z) = (0i*(O, . . . ,0;(O)- Then for ah t e T§°, 

i^t" I mi)-t-'ff{i\t)di < j^ti I {m)-t-'ffmdi. (7) 

.•1 Jc .-1 Jc 



i=l i=l 



where C = {l\li > ■ ■ ■ > /p > 0}. The right side of ([7j), the risk of is always finite 

( (fT8|) of Lemma [1]). Together with this finiteness, ([7]) leads to the inequality 

1=1 i=l 

where 

T,(t)=T,(ti,...,g=t^ / {Ul)-m)ff{mdl. (9) 



We also denote 

Tim(a; h) = Ti{ a, . ,a ,b, . ,b). (10) 

m p—m 

We will show that (IHl) implies Yl^=i ^«(*) = ^^d hence 0(Z) is almost surely equal to 
0*(Z) on C We integrate the both sides of ([H]) w.r.t. the measure t~^dt = (11^=1 t~^)dti . . . dtp 
over (0 < a < 6 < oo). Then we have 

5] / m<2j2 / / mii)-t-')tumt''dtdi. (n) 

The interchange of the integrals is guaranteed by (fT9|) in Lemma [TJ 



5 



Let 

= - + 1) mi) f nm f"-' fiM - f mmt-'^-^H^dt. (12) 

Then each integral of the right-hand side of f|TT]) is decomposed as follows; 

/ (^;(o-^.(z)) / 



where 



i?,(a, 6) = -i^ + 1) ' - Ul))G{l)H,{l; a, 6)rfZ, 



[ d,F{t\l) tidt + [ (^ + 1) F{t\l) u dt 

Ii{a, b) is bounded by the integral Ii{a, b) defined as 

/,(a, b) = K{^ + iy' - Ml)\G{l) 

d^F{t\l)r'^-H1dt+ [ (- + 1) F{t\l)r/''-Hidt 



dl. 



dl. (13) 



Lemma H] says that there exist constants Cim (^ = 1, . . . ,p, m = p — l,p) which are 
independent of a, 6 and satisfy 

p 

ii{a,b) < CimTlJ^{a;b). 

m=p—l 



Consequently with c = 2 maxj^m Qm we have 
P r p p 



J2 [ m)r'dt<cJ2 E Tll'{a;b)+2j2Ma,b) 

J 

i=l a i=l m=p— 1 1=1 

If we substitute and r (r > 1) respectively into a and b in dHj), we have 
j2 I m)t-Ut<cj2 E Tli\r-\r)+2j2Mr-\r), 

i=l -^"^(r) i=i m=p-l 1=1 



(14) 



(15) 




Figure 1: Relations among lemmas and Theorem 1 

where T(r) = T^„i. By Lemma fTOl there exists a constant M such that 
Tim{r~^] r) <M, 1 < V« < p, < Vm < p, Vr > 1. 

Since 

lim Ri{r-\r) = (16) 

r— >oo 

by Lemma [TTl the continuity of Ri{r~^,r) implies that Ri{r~^,r) is also bounded on the 
region r > 1 for each i. Therefore the left-hand side of f|T5l) is bounded and hence the 
increasing sequence limr^oo j'^r) -^i(^) converges for each i. This means 

/ Ti{t)t-^ dt <oo. (17) 
ixg° 

By Lemma fT2l the inequalities ( ITSl) . (|T6l) and (|T71) imply 

Ti(t) = 0, a.e. inl^, l<i< p. 

Hence ipi{l) = on C, 1 < i < p. (End of outline.) 

The following lemmas (see Figure 1 for their relation to Theorem 1) constitute a full 
proof of Theorem [H 

In the following Et[-] denotes the expected value w.r.t. the distribution of the eigen- 
values of Wishart matrix S with u degrees of freedom and the population eigenvalues 
(t^^, . . . , tp^)- We often use the inequality (x — y)^ < 2x'^ + 2y'^, x,y E M., to bound an 
integral from above. diag( ) denotes a diagonal matrix with diagonal elements 

Oi, . . . , Op. If Ai, . . . , Ak are square matrices of appropriate sizes diag(Ai, . . . , A^) denotes 
a block- diagonal matrix. 
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Lemma 1 

[{m)~tr^?fimdi<^ (18) 

i=l -^^ 

ill I mi)-umm)-t^Mm)dirHidt<^ m 

i=l -^"^a -J ^ 

Proof. The integral in (fT8|) can be written as ^t[Er=i(^**(0 - ^i^fif]- ^ suffices to 
prove Et[[ilJ*{l)f] < oo for each i. From ^ we have + 2)'ijj*{l) < YIj=ih = tr-S" and 
Et[{iYSf] < oo. 

Now we will prove f|T9|) . The integral can be written as 

Since the closure of is a compact region and the integrand is continuous in t on the 
closure of T^, it suffices to prove Y3l=\ -^t[lV'j*(0 ~ ~ < By Cauchy- 

Schwartz inequality, the following relationship holds. 

i=l 

<i2E,Mii)-t7') 

+ E i^^Mii) - t-r]V^' {E,[Mi) - 

i=l 

The last inequality holds since dominates I 
Lemma 2 Let a > 0, (3 > and b > a > 0. Then 

j-b / N 

where [x\ is the largest integer that is not larger than x and {a)j = a{a — 1) ■ ■ ■ (a — j + 1) 
is the falling factorial. 



1=1 i=l 

p_ 

i=l 

P 
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Proof. Note that {a)j > for < j < [a + 1] and {a)[a+2] < 0. By integration by parts 
I exp{—xP)x'^dx = ^(a°exp(— a/3) — 6"exp(— 6/3)) + ^ / exp{—x(3)x°'~^dx 

J a P P J a 



["+1] , N 

E exp(-a/3) - exp(-6/3)) 



j=0 ^ 



+ 4)1^ rexp(-x/?)x"-["+2]^^ 



["+1] / X 



Lemma 3 Let ), = {£c|(0 <)a < xi < ■ ■ ■ < Xp < b{< 00)}, Oi > 0, 

Pi > 0, i = 1, . . . ,p. Then 



J a;" exp ^— /?iXj^ (ia; 



is bounded by a linear combination of finite terms each of which has the form 

aELi(°»+7.+i) exp [ -a^ A / , 



with some integer vector 7 = (71, . . . ,7p). T/ie coefficients of the linear combination are 
positive and independent of a, b, (3i{i = 1, . . . ,p). 

Proof. By enlarging the region of integral to the direct product [a, &]^, we have 

exp(— /3iXi)(ixi X ■ ■ ■ X / Xp^ exp{—PpXp)dxp. 



/ exp — 22 Pi^i ] dx < x"^ 

Jx \ J Ja 



Applying Lemma [T] to each term on the right-hand side, we obtain the lemma. 
Lemma 4 Let Ii{a, b) be defined as in (|T3l) . Then for 1 < i < p, 



p 

ii{a,b) < CimTlj^{a;b) 

m=p—l 

with some constants Cim {m = p — l,p) which are independent of a, b. 
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Proof. 



X 



n^f \Fmt 



ti—ti + 1 



dL 



(20) 



where 



t'i (^1) • • • ; ti — ly • • • ; tp) y 

(!20|) is bounded by In + /j2, where 



K 



il 



'j2 



Z//2 + 

K 



Z//2 + 



t^i^'dt,di 



ti—ti + l 



&'di,dl. 



— ^i— 1 



First we prove the lemma for the case i ^ Let yg = Yl'j=ih^'sj- Then the inner 
integrals of In and Jj2 are rewritten respectively as 



n tf-'tr+iexpl-, 



^ ( Yl ^^y^ + ^^+1^^^ + ) ) dtidfi{H) (21) 

/ / n ^r^' ^-1 I E + + ^-i) ) ) di'MH). (22) 

If we use Lemma El the inner integrals of fl?I]) and fl2^ are seen to be bounded by linear 
combinations (whose coefficients are nonnegative and independent of a, b, ys{s = 1, . . . ,p)) 
of such terms as 



(a,+7,+l) 



JJyJ^expl-^ J2 Vsj 



where for (12T1) . 



and for ([22D, 



as = 1^/2-1, ys = ys, if I < s i,i + 1 <p, 
as = ys = yi + Vi+i, if s = i + 1. 



as = z//2-l, ys = ys, if I < s ^ i,i - 1 < p, 

as = ys = yi + yi-i, if s = i - i. 



(23) 



(24) 
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Consequently In and Ii2 are bounded by linear combinations of finite terms each of which 
has the form 

Ka^r<.^^<,(^s+,.+i) jj^^ ^ |^;(Z) -7A.(0|G(Z)neexp j d^x{H)dl, 



where a^, ^s, 1 < s 7^ i < p are given by (!23|) (for Jji) or (!24|) (for /j2), respectively. Besides 
the coefficients in the linear combination are nonnegative and independent of a, h. 
By Cauchy-Schwartz inequality, (|25ll is bounded by 



where 



[ [ (^*(Z)-V;.(Z))2G(/)exp(-^ y,U/i(^)rfZ, 

First consider the case for In. From (1251) . ([2]) and ([2]), we notice that 

A = T,p(a;6). (26) 
For the calculation of B, let X = (xij) ~ Wp(z/, S). We easily notice that 



B = E 



l[xl]^\^ = a~'l. 



where /„ is the p x p identity matrix and 



Xss ii 1 < s ^ i, i + 1 < p, 

Xss + Xs-is-i ifs = i + l. 



Therefore, with some constant K (independent of a,b), 

B = ka'^^^<-^^<p^\ (27) 
From ([23]), ([26D, and ([27D, it follows that ([25]) is equal to 

Now we consider the case for /j2. Similarly to the case In, 

A = T,p{a;b) (28) 
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and B = E [H,. x^J- |S = a-Hp] , where X = {xij) ~ Wp{u, S) and 



Xss if 1 < s 7^ 2, i — 1 < p, 

Xss + s+1 if S = Z — 1, 



hence 

B = Ka'^^^<^^^<p''\ (29) 
with some constant K (independent of a,b). From flM|) . fl28l) . and fl29|) . 

([25]) = iri/2a{^-»^-.^^+i>-''P/2-i^V2(^.^) ^ iri/2T//'(a;6). 

Finally we consider the case where i = 1 or p. Since the both cases are quite similar 
in the process of the proof, we only state a proof for the case i = p. li i = p the 
above argument for Ii2{= Ip2) still holds as it is and we only have to modify the part for 
hi{= Ipi)- The inner integral of Ipi equals 

[In ^f'' e^pf-^ E tsysY^'^'exp(-lyp)dt^d^^iH). (30) 
Jo(p) i<^.<p„i V ^ i<s<p-i / V ^ / 

By Lemma [3l the inner integral of fl30|) is bounded by a linear combination (whose co- 
efficients are nonnegative and independent of a,b,ys{s = l,...,p)) of such terms (the 
number of terms are finite) as 



a 



-i)5".+7.+i Yl y]^ exp (^-1 (^a 5^ + byp^ ^ 



and 



u/2-1 ifl<s<p-l, 
«s = < (31) 
I Z//2 II s = p. 

Consequently Ipi is bounded by a linear combination of finite terms such as 

Jc Jo{p) V,=i / 

X exp(^-^(^a^y, + 6|/p^^rf/i(if)rf/. (32) 

By Cauchy-Schwartz inequality, (15^ is bounded by 

a^YfsZl as+7«+l}-i'(p-l)/2^ap+7p-i'/2^1/2 ^1/2^ ^^33^) 
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where 



Similarly as before it turns out that 

A = Tpp_i{a; b), B = i^a-^^-' ^^r^^, 
where is a constant independent of a, b. Consequently (1331) equals 

i^i/2^{E?i^.+i}-KP-i)/2^a.-./2yV^2_^^^.^^ ^ ^1/2^;/^^ (a; 6). 

Lemma 5 For 1 < i < p, 

[ mi)-Ul)m) I F{t\l)r/'-\,dtdl<4Tl/\a,...,a) 
Jc J%i 

with some constant c- which is independent of a, b. 

Proof. Putting ys = Yl^j=i hh%^ see that the integral of the lemma equals 



■■{l)-Ul)\G{l) [ [ r/'-H,exp [-if^tsys] dtdi^iH)dl. 



By Lemma [3], the most inner integral is bounded by a linear combination (its coefficients 
are independent of a, b) of the terms whose forms are 



^T^Ef=i(ft+7.+-/2)g^p 



where 



p 

a 



u \ ^ 

2 2.%- 

i=i 



1, if j = i. 



By Cauchy-Schwartz inequality, 

aE^=i(/3i+7i+!-/2) 




CJOip) 
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is bounded by A^^'^B^^'^, where 

A = f _ ^^(/))2^(f) f (-- tr HLH'alp] dfi{H)dl, 

Jc Jo{p) V 2 J 

with L = diag(/i, . . . ,lp), while 

S = a^P/2+2E,^=i7. /" f G{l)y^^exp(--tiHLH'aIp]dfi{H)dl. 

We notice from and that A = (l/i\:)Ti(a, . . . , a). Let X = (xi^) ~ Wp(z/, S), then 

which is independent of a, 6. I 

At this point we need prehminaries about a partition before stating the next two 
lemmas. We partition (1, . . . ,p) into k blocks; 







= K-^E 




S = Ip 


\-j=l 











1st block (mo + 1, • • • , mi), 
2nd block (mi + 1, . . . , m2), 

fcth block {jrik-i + 1, ■ ■ ■ , m^), 



(34) 



where 



mo = < mi < 1712 < ■ ■ ■ < iTLk = p. 
Let [i], i = 1, . . . ,p, denote the number of the block containing i, i.e., 

[i] = s, if and only if m^-i + 1 < i < m^. 

(s), t = 1, . . . , k means the group of all the elements which belong to the sth block, i.e., 

i G (s), if and only if m^-i + 1 <i < m^. 

We also use the notation rhg = mg — ms-i, s = 1, . . . , /c, for the block sizes. 

Lemma El and Lemma [7] are just needed to prove Lemma [HI However these lemmas are 
useful in themselves since they give the asymptotic distribution of multivariate exponential 
type distributions under the block-wise dispersion of population eigenvalues. 

Lemma 6 Let eachpxp orthogonal matrix H = (hij) be partitioned as ( l55l) . There exist 
positive numbers 6i and S2{< 1) which are independent of H such that every orthogonal 
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matrix H has a series of pair {is,js), s = l,...,a; which satisfy the following three 
conditions. 

1- 1 < isjs < P and [is] > [js]- 

2. /i?. >5i. 

Isjs — ^ 

3. If i {1 <i <p) is not contained in \^ [m[j^]_i + 1, i^], then h^^ > 1 — 52, 

l<S<ijJ 

where [s, t] means the interval of integers from s to t. 

Note that the lemma includes the case that = 0, where the third condition V .^r., h^- > 
1 — ^2 for all i (1 < i < p) is the only condition to be satisfied. 

We give a proof of this lemma and Lemmas [TfHl below in Appendix. 

We still assume the partition ( IMI) for the next lemma. In addition, we introduce 
another condition and notation for the lemma. Let A*^") = diag{X^^\ . . . , X^^), n = 
1, 2, ... be the moving parameter matrix and we suppose that each X'f'\i = 1, . . . ,p, n = 
1,2,...) is decomposed as 

, (n) _ An) M p(n) ^ (n) ^ 

and 

limef^=e.(>0), z = l,...,p, (35) 

n—>oo 

hm a[jV«[? = 0, l<[j]<[^]<fc. (36) 

Us is the invariant probability measure on 0{ms). Ds,'Ss,ds {s = l,...k) are the 
submatrix or subvector of 

D = diag(rfi, ...,dp) H = diag(^i d= {di, ...,dp) 

respectively defined by the above-mentioned partition rule. Vs means the region given by 

{ds = idi)i(z(^s)\dms-i+l < ■ ■ ■ < drris}- 

Lemma 7 Suppose z^(> 0), at, i = 1, . . . ,p are given so that ai > (m[j]_i — i')/2 (1 < 
i < p). We also .suppose bij{> 0), 1 < [j] < [i] < k and Cjj(> 0), I < j < i < p,[i] = [j] 
are given. Let 

1/2 

(37) 



^*=i / \]>[j] 



As n ^ oo, the integral 

xr^^-^exp(^-^tTH'THA^'''^^dfi{H)dt (T = diag(ti, . . . , tp)) (38) 
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2ci, 

(i-ms-i){j-ms-i) 



converges to 

X exp( -itr/f;,i:),/f,,H, )rf/i,(/f,,)c?d, X TT /" x^^'' e^^(-]-x^]dx. 
^ ^ ^ HX.]-^" ^ ^ ^ 

iiTo 2S a constant which is independent of ai,bij,Cij, while 

K = ]^er-(P-"^w)/^ e. = a. - m[,]_i/2 + z//2 - 1 = 1, . . . 

Lemma 8 %(f) (1 < i,j < p) «s a bounded and scale-invariant function on C = {l\li > 
■■■>lp>0}- 

Lemma 9 Et[{litiY] is bounded intE T^. 

Proof. In the proof of Lemma 1 of Takemura and Sheena (2005), it is shown that 

P{tih >x\t)< P(x^(p_,+i) > x), Vx > 0,Vt G 1^, 
where xl{p-i-^i) is a chi-square random variable with u{p — i + l) degrees of freedom. Then 



POO POO 

Etiiktif] = 2 / xP{tik > x)dx < 2 / a;P(x^(p_i+i) > x)dx = E[x 
Jo Jo 



^ 1 

u{p-i+l)i 



Lemma 10 Ti{t) , i = 1, . . . ,p, are bounded in t E %\ 

Proof. First notice that 

p p p 



oo 



1=1 i=l i=l 

i=l i=l 

< AY,t'Mmi)-t^'f] = AY,Emm - 



i=l i=l 



The last inequahty holds since if) dominates t/?*. Therefore it suffices to show that 
Et\{'ip*{l)tiY] is bounded in t. From ([5]), we have 

From Lemma [H] and Lemma [9], the expectation of the right-hand side is bounded. i 



16 



Lemma 11 



lim Ri(a, b) = 

(a,fe)^(0,oo) 



Proof. Let Hi{l; a, b) be defined in (1121) . Using the monotone convergence theorem, we 
easily notice that Hi{l; a, b) converges to zero as (a, b) (0, oo). Clearly 



(^;(Z)-^,(Z))G(/)i/i(Z;a,6) <2|^*(Z)-^,(0|G(Z) / \d,F{t-l)\t^'^-Hldt. 
If the integral 

f \rdl)-Ail)\G(l) f \d,F{t-l)\t'''^-H1dtdl (39) 
is finite, then by the dominated convergence theorem 

lim Ri{a,b) = -K C^ + lY^ lim [ - ilJi{l))G{l)H^{l; a,b)dl 

= -K('^ + iy' [ mi)-m)Gil), ,lim HS;a,b)dl 

\Z / j£ (a, 6)^(0, oo) 

= 0. 

We will prove that fl39l) is finite. It suffices to show that the following integral is bounded 
in r > 1: 

/ \rAi)-Mi)\Gii) [ \d,F{t-i)\r'^-H1dtdi, 

where T(r) = {t|r^^ < ti < . . . < tp < r}. Since 
|9,F(t;Z) 1^/2-4^ dt 

= - f iutf") r\mt;l))tf-'dUdt, 

= {Hh(= r-i) < «i < ■ • • < < < • ■ ■ < ip < Vi(= I-)} 
the following equation holds. 

\i,l{i)-Ui)\G{i) [ \d,F{t;i)\r/^-HUtdi 

c Jl{r) 

ti=ti + l 

dt'- dl 

— — 1 



\m)-ui)\G{i) / r'^-' F{t\i)tf^' 

C Jl'ir) '- 



+ + [ \ijKl)-Ml)\G{l) I r'^-H,F{t\l)dtdl. (40) 



17 



The first integral on the right-hand side of poj) is bounded by 



mi)-i:m\G(i) 



n 



u/2-l 



l]L(r-\r) (seem 



mix 



u/2+l 



ti—ti + \ 



dt. 



ti — ^i — l 



dl 



(41) 



and the right-hand side is bounded in r > 1 by Lemma H] and Lemma [TOl Similarly by 
Lemma [5] and Lemma [TOl the second term on the right-hand side of ( l40l) is bounded in 
r > 1. I 

Lemma 12 The inequalities ( ITSj) . (fT6|) anc? (fT7|) imply 

Ti{t) = 0, a.e. m1^, l<i< p. 

Proof. We consider the terms on the right-hand side of (fT5|) . Fix i {1 < i < p) and 
m = p — 1 or p. Consider the following change of variables t — > {x,r) in each integration 
in (|T71) . where r and defined as 



r = 



Xq 



tv — s 



-s+l 



tp—s-\-l 

1 



s = 2,...,p— 1, ifm = p- 
p — 1, if m = p. 



(42) 



Conversely 



ti 
ti 



xi ■ --Xp-ir- 
xi ■ --Xp-ir- 



-\t. 



xi ■ ■■Xp-2r~ 
Xi- ■ ■ Xp-2r~ 



; tp—l 



X\T , tp 
X\T ^ tp 



r, ii m = p — 1, 
r~^, ii m = p. 



(43) 

We denote t expressed in terms of x and r in fH3|) by t{x, r; p—1) and t{x, r; p) respectively 
for the cases m = p — l,p. The domain of integral is shifted to 



< a;^ < 1, 
< X, < 1, 



s = 2, . . . ,p- 1, 
s = l,...,p- 1, 



< Xi < r^, if m = p - 
< r if m = p. 



(44) 



We can easily notice that Jacobian, J{t {x,r)) is given by 



,-p+i TTP"^ x^"-*^ 



-1 TTP-l 



if m = p — 1 , 



riLi a;? ^ if m = p, 



and that 



From 



p-2 rrp-1 ^s~p 

S 5 



^-1 



P-1 2,s-p 



if m = p — 1, 
if m = p. 



(45) 



(46) 



45|) and fl46l) , we have for m = p — 1 , p 



T,ft) t^^dt 



s=l 



Im{x,r)Ti{t{x,r;m))r drdx, 



(47) 
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where the indicator function Im{x, r) is given by 



I{xs <l,2<s<j9 — l)I{xi < r^), if m = p — 1, 
I(xs < 1, 1 < s < p — 1), if m = p — 1. 



For a while, we consider an inequahty with respect to Tj(t). We decompose Tj(t) as 



where 



/ 1 tsrls2hl^s2 ] dfi{H)dl 

J0(p) \ ^s,=ls,=l J 



X / exp I -- > ^ > ^ ts,lsX,s, 1 (4J 
For the two points 

t^"^^ = t{x^^\r; m) = t{x^i \ . . . , x^\r; m), t*-^-* = t{x^'^\r; m) = t{xf\ . . . , x^\r; m) 

defined by (H3l) with x^^\ x^"^^ such that x^"'^'* < x^p {j = 1, . . . ,p — 1), we have the following 
inequality 

'') 

Notice that T/^^(t(2))/Tf is independent of r, since it has the form YYj=i{^f'^ 
with some constant a/s. 

Let J\f = {x\c < Xj < 1, j = 1, . . . ,p — 1} with some constant < c < 1. If we apply 
the inequality (B9|) to the two points 



= t(a;,r;m), x G AT, t^^'^ = t{l,r;m), 1 = 

p-i 

we have 

Ti{t{l,r;m))<R,m{x)T,{t{x,r;m)), Vx e Af, (50) 

where 

. ._ Tf^(t(l,r;m)) 
Tl^>{t{x,r;m)) 
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Now we evaluate integral fj47|) using the inequality flSOjl . Since Tj(t(l, r; m)) = Tjm(?" ^ 
for any x G A/", 

/OO /"OO 
Ti^(r-';r)r-Mr < J Rimix)Ti{t{x,r;m))r~^dr. (51) 

Notice that if as G A/", then 

J^(l,r) = J(r > 1) < J(r > x}/') = /„(a^,r) 
and the compactness of J\f implies that there exists some c*(> 0) such that 

Rimix) < C*. 

Combined with (I3T|) . this means that for Vx G A/", 

/OO roo 
%^{r-^] r)r-'dr < c* J I^x, r)%{t{x, r; m))r-^dr. (52) 

Suppose that there exist 5{> 1) and e(> 0) such that 

T^„^(r-V) > e forVr>5, 

then 

/•OO /'OO 

/ Tim{r~^;r)r~^dr > e / r^^rfr = oo, 

J(5 Js 

which implies that the integral on the right-hand side of fl52|) also diverges. This fact to- 
gether with P71) implies J^oo Ti{t)t~^dt = oo, which is a contradictions to ([T7|) . Therefore 
we can conclude that for any 5(> 1) and e(> 0), there exists r such that r > 6 and 

Tim{r~^;r) < e. 

This enables us to construct a series rj{j = 1, 2, . . .) such that oo and 

T,„(r-Vj)^0 (53) 
as j —>■ oo. This folds for any i {1 < i < p) and m {m = p — l,p). From ([TB|) . we have 

lim Ri{rT\ r,) = 0, 1 < Vi < p. (54) 

J— >oo 



It follows from ([l5]), (153]) and ([5lD that 

P r 

V / T,{t)t-^dt = Q. 

^=l 

Therefore Tj(t) = 0, a.e. in T[f , 1 < ^ < p. 
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5 Appendix 



Here we give proofs of Lemmas El [7] and [HI Correspondingly to the partition stated before 
Lemmas El we make the following partition of a. p x p matrix A = {aij); 



/All ■■■ Aik \ 



A 



Ast '■ rris X nit matrix, 1 < s,t < A;. (55) 



y Aki ■ ■ ■ Akk j 

For a vector a = (ai, . . . , Op), the corresponding partition is given by (ai, . . . , a^). 

Proof of Lemma [6l We use the notation . . . ,t) (s < t) as the principle submatrix 
that consists of the blocks Hij, s < i,j < t. Namely H{s,...,t) consists of all the 
elements hij such that s < [i], [j] < t. From now on if we refer to a "submatrix", it only 
means a principle submatrix that consists of the blocks. 

Choose a small enough positive number 6o- We define the term "separable" with 6o- 
Consider a submatrix H{si, . . . , Sp). If for some i {I < i < p), the squared sum of 
blockwise-off-diagonal elements J2si<[j]<s, Si+i<[{\<sp^'ij smaller than 6o, we call this 
matrix "separable" (into H{si, . . . , Si) and i/(sj+i, . . . , Sp)). If we make a repetitive 
separation, starting with H itself, finally we have a series of submatrices (not necessarily 
unique) 

H{l,...,si), H{si + 1, . . . , S2), ■ ■ ■ , i3"(sK_i + 1, . . . , s«;), 

{1 < Si < ■ ■ ■ < = k), each of which is unseparable. Though these matrices are not 
necessarily orthogonal, if the lemma holds for each of them, obviously it also holds for H 
itself. We easily notice that there exists a positive constant c (independent of H) such 
that 

^i<[i]<si > 1 — c5o for any i such that !<[«]< Si, 
Ei<[ii<si > 1 - c5o for any j such that 1 < [j] < Si, 



I]s,_i+i<[j]<s, hfj > 1 - c5o for any i such that + I < \i] < s^,, 
.Es„_i+i<[i]<s„ ^1 > 1 - c5o for any j such that s«_i + 1 < [j] < s^. 

Therefore we only have to prove the lemma for H under the condition that H is not 
necessarily orthogonal but unseparable and satisfies the conditions 

J2 hi > 1 - c5o, 1 < Vz < p, (56) 

i<i<P 

J2 hi > 1 - c6o, 1 < Vj < p. 

l<i<p 
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First consider the case k = 1, namely ii" is a single block matrix. If we put 62 = {c+ l)6o, 
then (l56l) imphes 

E^'.= E hl>l-S2, l<W^<p. 
je[j] i<i<p 

The lemma holds as the case u = 0. 

Now we suppose k > 2, where H consists of multiple blocks. First since H is unsep- 
arable, we have 

[i]=fc,l<[i]<A;-l 

which means there exists some z's (g (k)) such that 

E ^ Sorhk^. 

i<[i]<fc-i 

Put the largest i as among z's that satisfy (jS]). Furthermore guarantees the existence 
of j (1 < [j] <k-l) such that 

Put this j as ji. The way ii is chosen implies 



hi < Som^^, if ii < i < p, 

i<[i]<fc-i 



which means if ii < i < p, 



[j]=fc i<b1<fe i<[i]<fc-i 



> E ~ ^O^k ^ = E ~ ^O^k 
l<[j]<k l<j<p 

> 1 — c5o — ^o'lTT'k^ (because of (l56l)). (57) 
We proceed to the second step. From unseparability of H, 

E ^0. 

[ii]<[j]<fc,i<[i]<[ii]-i 

This means for some such as [ji] < [i] < k,l < [j] < [ji] — 1, 

> 6omyl^_^{p - my^yi)-\ (58) 

Choose the smallest j (and, if necessary, the largest i) among (i, j)'s that satisfy fl58l) and 
put these i,j to be 12, 32 respectively. 

Repeat the "second step" until [js] reaches 1. (Note that [js] is strictly decreasing as 
the step is repeated). Finally we have a series of {is,js), s = 1, . . . ,u, where [ju] = 1, 
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and [js] < [is+i] {s = 1, . . . ,uj — I). Obviously {is,js) satisfies the first condition of the 
lemma and the condition 



hi^j^ > ai, ai = 5o"^fe 

hlj^ > a„ a, = 6omyl^^^_^{p - my^_,]_i)~\ s = 2,...,u. 

Note Qs > SqP^'^, 1 <Ws < uj. If we define 6i as 6i = 6op~'^, then the second condition is 
satisfied. 

Finally we consider the third condition of the lemma. Notice that [jg] < [ig+i] implies 
my^]_i + 1 < is+i, hence 

[1, p]/ 1) [m[j^]_i + 1, is] = [{max is) + 1, p] C [ii + 1, p]. 

1<S<UJ 

Therefore if i ^ [Ji<s<uj[''^[js]-i + "^s], then ii + I < i < p. From (1571) . for such i, 

[j]=k 



since 62 = {c+ l)So > {c + irtf, 



0- 



Proof of Lemma [T], In a small neighborhood, an orthogonal matrix H has its 

strictly (left-)lower part {hij)i^j as its coordinate function; H has one-to-one correspon- 
dence to {hij)iy,j, and {hij)i<j is a C°° function of {hij)iyj. Since 0{p) is compact, 
we have a finite coordinate neighborhoods, {0^'^\ (f)r) , t = 1,...T for 0{p) such that 
4>t{H) = u = {uij)iyj, Uij = hij{i > j) for H G O^'^^ or conversely 

[hjj{u), iil<t<3<p, 

where hjj{u) is C°° function on Ur = (priO^'^^). 

Let Jt{u) denote the Radon-Nikodym derivative of /i with respect to the Rp^P'^^^- 
dimensional Lebesgue measure, i.e. JT-{u)du = d^{H). Actually 

d^I{H) = co/\{h,ydh,, 

i<j 

where hi {i = 1, . . . ,p) is the ith column of H and Cq is a constant. If we build in fl59l) 
and the fact 

duij if 1 < j < i < P, 



dh 



T.s>t o^,dust ifl<i<3 <p, 
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into the above wedge product, Jr{u)du is obtained. If we use the partition of unity Lr{H) 
subordinate to O^'^^ (r = 1, . . . , T), the integral fl38|) is rewritten as 

X JJ u^' JJt^/'~'exp('-itr(/fW(n))'T/f(^)(w)A("Adwdt, (60) 

where H^'^\u) = (hij) is given by (!59|) . 

Consider further change of variables {t,u) {d,q), where d = {di, . . . ,dp) and 
Q = {(lij)i>j, given by 

di = tja^.j 

^ juij if i > j, [i] = [j], 

The Jacobian is given by 

J{{t,u) ^ (d, q-)) 
= J{t d)J{u q) 

= ri(«[?)" n <^"'wr""^i4\y'" 

i=l i=l j=l [i]>[i] 

Notice that 

n(M,<"')°" = n(*«.'"')'" = n<- n(«.'"')°- (62) 

i=l i=l i=l i=l 

\{tr-' = \{d'J'-'\{{a^y''^' (63) 

j=l 4 = 1 j=l 

n = n (64) 

W>[j] W>[i] 
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From ([61]), ([62D, dM]) and dMD, the integral §0^ equals 

P ^ 



1=1 [i\>[j] i>j,li\ = lj] 

^ s=l i.ie(s).i>i i.ie(s).i<i ' 



k 

i I ^ 

X exp 



i,j&{s),i>j i,j&{s),i<j 



+ E4+E('4''M)<H"[,>[?) 

W>[i] W<[i] 



dqdd, (65) 



where i?+ is the positive part of R, Ci = ai — (p — m[j])/2, Cj = aj — m[j]_i/2 + z//2 — 1, (i = 
1, . . . ,p) and Ijy(„){d) is the indicator function of the region 

V^-^ = {d\d^{a^^r<...<d,{a\;lr}. 

The notation u in the integrand is an abbreviation of n(q, d, ck*^")) (^*^") and cx^^^ 

respectively means ^"^"^ = (^j"\ . . . , ^p"^) and a*^"^ = (q;i"\ . . . , a^"'')) which is specifically 
given by 

{Qij ii i > J, [i] = [j], 

In order to evaluate ( l65i) , we use Lemma El By the lemma, every orthogonal matrix 
H has a set of pairs {is,js){s = 1, . . . ,uj) that satisfy the conditions of Lemma [61 Define 
Tij{H), 1 < j < i < p as an indicator function as follows; 



1 if (iJ) = (isjs), 1 < 3s < a;, 
otherwise. 



Then every H has an index of {Tij{H))i<j^i<p. Since the existence of {is,js){s = 1, . . . , cu) 
may not be unique, H can have more than one index number. However if we put a prefer- 
ence order among all possible (2^^*'"^)/^) index numbers, the index is uniquely determined. 
By this index, we can naturally partition 0{p) into the subsets O^'^^ (f = 1, . . . , 2^^^"^^/^). 
Let the corresponding partition of unity be denoted by rf(i3"). Now f[^S]) is expressed as 

p T p(p-l)/2 

n(«.'"V'E E 

j = l T=l f=l 
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where Irf is given by 



p 



X 



n< n ^ n ^ 

i=l [i]>[j] i>j,[i] = [j] 

s=l i.ie(s).i>j i.ie(s).i<i 



i,j&{s),i>j i,j&{s),i<j 

, f(n)/ ,(n) / (n)^ 



+ E4+E("S;'("))''4"'("S'/"S') 

H>[i] [i]<b] 



dqdd. 



Now we focus on I-rf. Take large enough n. Suppose (d, q) satisfies 

Then > imphes H^'^\u) has a sequence {is,js), s = 1, . . . ,uj that satisfy the conditions 
of LemmaEl First suppose i G Ui<s<(^["^bs]-i + l''^^] (^^-Y -^(^))) then for some s, 



681) implies 
691) is equivalent to 



(68) 
(69) 

(70) 
(71) 

(72) 



Moreover the fact /^icn) > implies 
From (COl), (IZID and (I72D, if i G /(f), then 

(n) (n) 

"[i.l "[is] [ii]>[i2] 

where in the third inequality we used the fact there exists a positive number ^ such that 
C < ^1"-' for alH (1 < i < p) and all large enough n. Consequently 

Ifid,q) = l (73) 

under the condition (l67j) . where /f (d, q) is the indicator function of the region 



\{d,q) d,<6^'r' E V^G/(f)|. 



[ii]>b2] 
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On the other hand, if i ^ -^(t); the condition 3 of Lemma [6] guarantees 



which means 



s=i \,je{s),i>j i,je(s),i<j 

= E*( E 4+ E 



(74) 



M-j- = < 1. 



We also notice that l^- > imphes that if i > j, [i] = [j], then 
Therefore under the condition (1671) 

^QiQd) = 1, 



(75) 



where iQ^Qd) is the indicator function of Qd = {<lij)i>j,[i]=[j] with respect to the region 
< 1, 1 < j < < = [j]}- From ( 1731) and ( 175|) . the following relations hold. 

Since Jt{u) is bounded on a compact set, the integrand of Irf is bounded by 

j=i [j]>j] «>j,[«]=[i] 



X exp 



H«E( E 4*+ E (".'I'M)'*) + E ««} 

s=l \,je{s),i>j i,j&{s),i<j ^ [i]>[j] 



(77) 



with some constant c. 

From (17il) . ( 1761) and (177j) . we notice that the following function f{d,q) dominate the 
integrand of Irf] 

7{d,q) = ch{d,q)lQ{q,)fldt' H qf^ H 4"^' 

i=l [i]>[i] i>J,[j] = [i] 

X exp|-^^(l - ^2) ''''P{"^ 5Z 



H>b1 
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We have 



R+P jRp(p-i)/^ 



f{d, q) dq dd 



i>j,[i] = b] 



i0{f) 



X 



■^^''^ ie/{f) J H>i] [i]>[j] > 



(75 



where d^^) = d^^) = {di)i^i(f), = #{1 < ^ < p|z ^ /(f)}, = #{1 < ^ < 

p|i G /(t")}. Since Cj > —1, 



i0(f) 

is finite (say Mq) and independent of q, while 
II df hid,q)dS'^= II 



ie/{f) ie/(f) ^ [i\>lj] ^ 

ie7(f) Jo \ J 



< 00 



It turns out that ( ffSl) equals 



cMo JJ Mi (5i^)-E.e/(f)(^'+i) 

jG/(f) 

n4'"(E4 



X 



^''^ [i]>[i] ^W>[j] 



^ W>[i] ^ 



X 



/ n lQ{Qd)dqd, 

jRPl . 



i>j,[i\ = [j] 



where qo = {qij)[i\>\j] and pi = Y!l=i^s{ins - l)/2, p2 = Y.[i]>[j] mmj. This integral is 
obviously finite since % > ([z] > [j]), + 1 > (1 < z < p), Cij > (z > j, [z] = [j]). 
The finiteness of (1751) guarantees the use of dominated convergence theorem. Therefore 
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as n oo, Irf converges to 



R+P jRPip-^)/^ 
k p 



X 



X exp 
where 



U'^jA) n<- n n 

s=l i=l [{\>[j] i>j,[-i] = [j] 

^ s=l i,j£{s),i>j i.iG(s).i<i ^ li]>\i] ^ 



dqdd, (79) 



u* = lim u(q,d,^("),a(")) 

n— >oo 



= lim {ud,Uo) (u = {ud,Uo), Ud = {uij)i>j,[i]=[j], Uo = {uij)[i]>[j]) 

n — s-oo 

= ( lim Ud, lim Uo) = {qd, 0) (see dSS}) 

n— >oo 

and Ix>Xds) is the indicator function of the region 

T^s = {dg = idi)i(=(^s)\dms-i+l < ■ ■ ■ < dm,}- 

If we change the notation as qd Ud, (ff9l) equals 

Lr{H^^\{ud,0)))l,{H^^\{ud,0)))M{ud,0)) 




k p 



X 



n^^=(^^) n n 

s=l i=l [i\>[j] i>j,[{\ = [j] 

i{E( E (H'"(K.o)))j,d.e 



X exp 



i,j<^{s),i>j 



[i]>b] 



dqo dUd dd. 



^0) 



^{p) = {H\Hss (1 < s < fc) are all orthogonal matrices} is the subgroup of 0{p). 
H G 0{p) if and only if all the off-diagonal blocks {Hgt-, s ^ t) are zero. Therefore we 
can identify 0{p) with the product group (9 (mi) x ■ ■ ■ x 0{ffik)- Notice that 



/f = ifM((n,,0)) 

is in 0{p). 

Consider the following transformations of H G 0{p); 

H HH\ 



(81) 



^2) 
^3) 
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where H* = dia.g{Hl-^^, . . . , Hlj^) G 0{p). If we consider these transformations on 0{mi) x 
■ ■ ■ X 0{mk), they are equivalent respectively to 



Hgs ^ HggHg,,, s — 1, . . . , k, 

Hgs ^ -^ss-^ss^ s = 1, . . . , k. 



^5) 



The unique invariant probability measure on 0{mi) x ■ ■ ■ x 0{mk) with respect to the 
both transformations (IMl) and ( l85l) is 

/il X ■ ■ ■ X /ifc, 

where /i<j is the uniform probability measure on 0{ms) {s = 1, . . . ,k). 

Now we examine the measure on O^'^'* fl 0{p) given by Jr{{Uii,0))dUd through flSTl) . 
JT-(('iA(i, 0))(iitd is derived from Jr{u)du by imposing the condition ito = 0. We easily 
notice that under the condition Uo = 0, cq /\^^-{hi)' dhj{= Jr{u)du) equals 



s=l l<i<j<ms 

where hf'\i = 1, . . . m^) is the ith column of Hgs. This differential form gives the invariant 
measure on 0{p) with respect to the both transformations fl82|) and fl83|) . Therefore 
Jr{{ud,0))dud gives the invariant measure (w.r.t flH^ and fl85|) ) on 

{(/fn, . . . , /ffefc) G 0(mi) X ■ ■ ■ X 0(mfc)|diag(ffn, . . . , H^k) G O^^^}, 

which is equal to Kq/ii x ■ ■ ■ x /i^ with some constant Kq. (Note that Kq is independent 
of ai,bij,Cij,T.) Consequently flHOl) is equal to 

Ko / / ■■■ / / Lr{diag{Hu,...,Hkk))lf{dmg{Hii,...,Hkk)) 

Jr^_^_ Jo{fhk) Jo{mi) Jrp2 



X 



s=l j=l [i]>[j] s=l [i\ = [j]=s,i>j 

k \ / ^ \ k 



X expf-^'^ti H'^^DsHss^sj X expf-i ^ q-Adqo f[dfisiHss) dd. 
Adding up (l86l) over all r's and f 's, we have 



(86) 



, 2c,: 



X exp 
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Proof of Lemma [81. The property of scale-invariance is obvious from its definition. 
We will prove the boundedness. For the case i < j, it is obvious; 

njH) = Tij{l)^-f < < ^ 



Now we suppose that i > j. We notice 



[^ + 2)%(Z) = ||^, (87) 



where 



hil)= [ [ ihlul,){Uh)r/^-'exp(-ltiH'THL]dfi{H)dt, 
Ji"^ Jo{p) V 2 / 

JO(p) \ ^ J 

with T = diag(ti, . . . , tp), L = diag(/i, . . . ,lp). Furthermore we notice that /i(Z)//2(Z) is 
equal to /i(Z)//2(Z), where 

h{i)= [ [ ihltJ,){Uk)r/^-'exp(-ltiH'THL]dfi{H)dt, 
^•X5° Jo{p) V 2 / 

Ui)= I I itiliff/^-'exp(-l-tiH'THL]dfiiH)dt, 
Ji-^ Jo(p) \ ^ ) 

with It = lt/k(t = 1, . . . ,p)J = (!i{= 1), 12, . . .J Ip) and L = diag(/i(= 1), ■ ■ ■ , Ip)- 

We will prove that /1//2 is bounded on £ = {Z|l > I2 > ... > Ip > 0}. First let 
(^2, • • • , Ip) be parameterized as follows; 

rt = lt+i/h, t = l,...,p-l, 

equivalently 

i-l 

k = Y[rs, t = 2,...,p. (88) 



s=l 



C is equivalent to 7?. = {r = (ri, . . . , rp_i)|l > > (t = 1, . . . ,p — 1)}. It suffices to show 
that /i(Z(r))/J2(Z(r)) is bounded on 71, where Z(r) is given by (IHHj) . It is easily proved 
that /i(Z(r))/J2(Z(r)) is continuous on 71. If we can expand Ji(Z(r))//2(Z(r)) continuously 
over 7Z = {r = (ri, . . . , rp_i)|l > > (t = 1, . . . , J9 — 1)}, then the expanded function 
is continuous on a compact region 7Z, hence is bounded. Therefore we only have to show 
that for an arbitrary sequence r^^^ = {r^\ . . . , r^^) and r G 7Z such that 

lim rj"^ =rt, 1 < t < p. (89) 

n— >oo 
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/i(/(r("^))//2(Z(r("^)) converges to a certain value which depends only on r. 

Now choose an arbitrary sequence r^"-* = {r^i \ . . . that satisfies (1891) . Let /c — 1 

denote the number of r^'s in (!89|) that are equal to zero. If > 0, 1 < Vt < p — 1, then 
the continuity of Ji(Z(r))//2(Z(r)) on 7?. guarantees 

hm /i(Z>(")))//2(/>("))) =/i(Z>))//2(Z>)). 

n^oo 

From now on we suppose that k > 2. We define ms{s = 0, . . . ,k) so that 

frui = 0, i = l,...,fc 1, 

"^o(= 0) < mi < . . . < mfc(= p). (90) 

Based on the partition (MI]) , we apply Lemma[7]to /i(Z(r(")))//2(Z(r("))). Let Z = {slr^ = 
0, l<s<p — 1} = {s\s = rrit, 1 < 3t < A; — 1}. Let lt(r^"'^) be decomposed as 



s=l 

where 



l<s<t-l,se2 l<s<t-l,s^Z 

Notice that if [t] = u, 

u-l 

n 



^j."^ and Oj^"'' satisfy the conditions (l35ll and ( |36ll respectively; 



^l<s<t-l.s^2 ^ 



lim 

n 



hm = lim J] =0- 1 < [^2] < [h] < k. 

S=[t2] 

/i(Z(r(")))//2(Z(rW)) equals I^/I^, where 

and is given by (1371) . Now we can apply Lemma [71 First if [i] > [j], then as n — 00, 
II converges to 



X 

[tl]>[<2j 
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where 

_ I J~j'^-(P-™[t])/2 j^l-(p-m[i])/2 

1 - m[i]_i/2 + Z//2 - 1 if t = z, 
-m[t]_i/2 + Z//2 - 1 if t 7^ z, 

1 if ti = z, t2 = j, 



^tit2 - ^ (91) 
otherwise. 

If [i] = [j] , i > j, we rewrite /* as 
ft-(r("))A(r('^)))(K("))-i 

The fact ^ ^ 

lim M^^^l = liin ^ = |l 

and Lemma [7] imphes that converges to 

k 



Ko K n / / \{dT n (^-) 

X exp|^-itriJ^^i:)si?ssSs^d/is(J?ss)c?c?s X |^y" exp da;^ 



(ti-ms_i)(t2-ms-i) 



where 

-(p-r7i[t])/2\ 2-(p-m[,])/2 



K = ( ]-[^-(P-WV2 W2-(P-H)/2^ ^g2) 



^ ^ I 2 - m[,]_i/2 + Z//2 - 1 if t = z, 
\-m[,]_i/2 + z//2-l ift^^, 

with &f^f2 as in (1911) . On the other hand, using Lemma [7] again, we notice that I2 converges 
to 

2 ; ) 



te(s) 

00 / 1 \ \ J2s>t fhsfrit 



X ( / exp 
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where K and are respectively given by fl92l) and fl93|) . Consequently we notice that in 
either case, 11/ 12 converges to a certain value which is dependent only on r. | 
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